Thermodynamics of a Fermi liquid in a magnetic field 
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We extend previous calculations of the non-analytic terms in the spin susceptibility Xs(T) and 
the specific heat C(T) to systems in a magnetic field. Without a field, Xf{T) and C(T)/T are linear 
in T in 2D, while in 3D, Xs(T) oc T 2 and C{T)/T oc T 2 logT. We show that in a magnetic field, 
the linear in T terms in 2D become scaling functions of ^lbH/T. We present explicit expressions for 
these functions and show that at high fields, hbH S> T, Xs(T, H) scales as \H\. We also show that 
in 3D, Xs{T, H) becomes non-analytic in a field and at high fields scales as H 2 log \H\. 
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Landau Fermi liquid theory provides the basis for 
our present understanding of correlated electronic sys- 
tems. The theory predicts that, in any Fermi liquid, 
the spin susceptibility Xs(T) and the specific heat co- 
efficient 7(T) = C(T)/T tend to a constant at T -> 
P, Q]. Later, Landau theory has been extended to in- 
clude the leading temperature dependence of Xs{T) and 
7(T) which turn out to be non-a naly tic in dimensions 
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zero-temperature terms, the thermal corrections come 
from fcrmions in the immediate vicinity of the Fermi sur- 
face. In 2D systems, both Xs(T) and j(T) are linear in 
T 0, @, 0, El and the coefficients are expressed in terms 
of charge and spin components of the scattering ampli- 
tude at the scattering angle 9 = ir Jl2j j l3 | . In 3D, Xs{T) 
is quadratic in T, i.e., is analytic |4j5U6|, while "f(T) is 
non-analytic and scales as T 2 log T 0, Q • 

In this communication, we extend previous works to 
systems in a magnetic field H. We consider S = 1/2 
charge-less fermions (like 3 He atoms) for which the mag- 
netic field adds spin-dependent Zeeman term ±/is_ff to 
the fermionic dispersion. We show that, in the pres- 
ence of a field, A Xs (T,H) = Xs(T,H) - x«(0,0) and 
Aj{T,H) = 7(T, H) - 7(0,0) become scaling func- 
tions of fj, B H/T: A Xs (T,H) = A Xs (T,0)f x (» B H/T), 
Aj(T,H) = Aj(T,0)f 7 {fj, B H/T). We present the ex- 
pressions for these functions to second order in the in- 
teraction potential U. For 2D systems, we show that at 
[i B H > T (but still, fi B H < E F ), Ax s (T,H) scales 
as \H\ and weakly depends on T. In the same field 
range, 87 (T, H) is still linear in T, but the prefactor 
is different from that at H = 0. For 3D systems, we 
show that A Xs (T,H) becomes non-analytic at a non- 
zero H. The non-analytic term in Ax s (T,H) scales as 
H 2 \og[max(n B H, T)/E F ]. The specific heat coefficient 
j(H,T) in a field still scales as T 2 logT, but, like in 2D, 
the prefactor changes between H = and [i B H ^> T. 

The analysis of the behavior of A Xs (T, H) and 
A"f(T, H) in a magnetic field may be useful for experi- 



mental verifications of the non-analytic behavior of ther- 
modynamic variables. It is more straightforward to an- 
alyze the dependence on the magnetic field rather than 
the dependence on the temperature. In particular, recent 
measurements of the temperature dependence of the spin 
susceptibility in Si inversion layers |l5| didn't yield con- 
clusive results on whether the T dependence of X s{T) 
is indeed linear, as some temperature dependence in- 
evitably comes from spins on the substrate. We propose 
to measure the field dependence of the spin susceptibility 
at a given T and use our scaling functions to fit the data. 

The point of departure for our calculations is the 
Luttiger-Ward expression for the thermodynamic poten- 
tial. To simplify the presentation, we assume that the in- 
teraction potential U (q) is independent on q. We restore 
the momentum dependence of U{q) in the final formulas. 
To second-order in U , the thermodynamic potential is 
given by 
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where $0 is the thermodynamic potential for free 
fermions, and Tl^(q,n n ,T) and JI ll (q,D, n ,T) are the 
particle-hole bubbles composed of fcrmions with spin up 
or spin down, respectively. 

Previous studies of the spin susceptibility and the 
specific heat in a zero magnetic field established that 
the non-analytic temperature behavior of A Xs (T) and 
Aj(T) originates from the non-analyticity of the polar- 
ij^atipii operator either near q = (Landau damping) 
lid] or near q = 2kp (a dynamic Kohn anomaly) 
The 2k F non-analyticity contributes 
to the spin susceptibility and the specific heat, while 
the q = non-analyticity only contributes to the non- 
analyticity in the specific heat. This can be easily un- 
derstood as the non-analytic term in the zero field spin 
susceptibility A Xs (T) describes a singular response to 
an infinitesimally small magnetic field. A magnetic field 
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FIG. 1: The diagram for A"l>2fc F , and the trick to compute it. 
As the non-analytic part of A$2fe F comes from small k and p, 
it can be re-expressed as a product of two bubbles lV^(q',u) 
with small momentum transfer q' = k — p. 



splits Fermi momentum kp into kL and fe*L The small q 
form of n T T (q, fi n , T) and II U (q, O n , T) is unaffected by 
this splitting, up to terms of order (fisH/Ep) 2 , hence 
the response to the infinitesimal field must be analytic in 
T. At the same time, singular 2kp contribution to $(T) 
at zero field originates from the fact that the two polar- 
ization operators in Eq. Q are non-analytic at the same 
q = 2k p . In a field the singularities in spin- up and spin- 
down polarization operators occur at different q = 2k F 
and q = 2k p. Accordingly, a magnetic field regularizes 
2kp non-analyticity in the thermodynamic potential, but 
for a price that the linear response to the field, i.e. the 
spin susceptibility A%(T, H = 0), becomes non-analytic. 

Our goal is to analyze the forms of the susceptibility 
and the specific heat at a finite H, i.e., beyond the linear 
response theory. We consider the fields for which /j,bH 
is comparable to T, but still ^bH <C Ep. For these 
fields, the non-analytic contribution to $ from small q are 
unaffected by the field . However the 2kp contribution is 
field dependent and evolves at fisH ~ T. 

The calculation of A$ = <f> — <£>o is somewhat tricky. 
In principle, all one has to do is to evaluate particle-hole 
bubbles for fermions with up and down spins at a finite 
T, substitute the results into Eq. (JIJ, integrate over mo- 
mentum q and sum over Matsubara frequencies f2„. In 
practice, however, this computation is easy to perform 
only for small q part as for q <C kp, the non-analytic 
part of the polarization bubble is associated with the 
Landau damping, which does not depend on T, apart 
from regular (T/Ep) 2 corrections. Accordingly, one can 
safely use the known analytical forms of II(<7, J7) at T = 0. 
For q near 2kp, non-analytic terms in (q,fl n ,T) and 
Il-'-^g, fi n , T) contain scaling functions of T/u>, which are 
only available in integral forms 0. This substantially 
complicates direct calculation of the 2k p term. There 
exists, however, a way to compute the 2kp term, which 
avoids dealing with the 2k p polarization bubbles at a 
finite T. This method explores the fact that only the 
non-analytic parts (q,Q n ,T) and U^(q,Q n ,T) for 
q near 2kp contribute the non-analyticity in the ther- 
modynamic potential. Earlier works have demonstrated 
that the 2kp non-analyticity in H(q, f2 n ,T) comes from 



fermions in the particle-hole bubble with momenta near 
±<f/2 This implies that, out of four fermions in 

the second order, two-bubble diagram for the thermo- 
dynamic potential in Fig. 1, two fermions with oppo- 
site spins have momenta near q/2, while the other two 
fermions have momenta near —q/2. Then the 2kp part 
of the A$ can be re- written as the integral over small k 
and small p of 
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or, equivalently, as 
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IlU(q',Sl n ,T) 



(3) 



where the integration is confined to small q' = k — p. 
In other words, the non-analytic 2kp contribution to the 
thermodynamic potential can be re-expressed in terms of 
the particle-hole bubble for fermions with opposite spins 
and a small momentum transfer. The non-analytic term 
in IT at small q' does not depend on temperature (apart 
from irrelevant corrections), hence H^(q',Sl n ,T) can be 
safely approximated by IV*(q', 57 n , 0). At the same time, 
the polarization bubble H^(q, Q n , 0) strongly depends on 
the magnetic field (contrary to (q,tt n ,0)), and this 
gives rise to the scaling dependence on hbH/T. 

Combining the q = and 2kp contributions, we obtain 
for the thermodynamic potential 
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where the integration involves only small q. 

We next proceed separately with 2D and 3D cases. In 
2D we have 



U^(q,n n ,T) = 
nTT(g,n n ,T) = 



I On I 



2n ^(n n - i Stf + (vpq) 2 
m I fi„ I 



2n ^n 2 n + (vpqf 



(5) 



where dots stand for analytic terms, expandable in pow- 
ers of f2„ or q 2 , and 5fx = /Xf — Hi = 2[1bH . Substituting 
Eq. (JSJ into Eq. Q and integrating explicitly over mo- 
mentum q we obtain 



. _ /m\2 U 2 T 2 , 
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(n„ - 2i f i B Hyn 2 n 



(6) 
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FIG. 2: Ax as a function of temperature at fixed mag- 
netic field (a) and as a function of magnetic field at fixed 
temperature (b). In the left panel, T and Ax/ fi%A are 
in units of {IbH (A is defined in the text). In the right 
panel /ibH and Ax/fJ.%A are in units of T. In these units, 
A X /AA = f x (VBH)). 



Differentiating Eq. JHJ with respect to H, we obtain 
<9A$ '■- 4 - 4 



AM 



dH 



V n n 2 + (2fi B H) 2 



n 3 k 2 F 



The sum over Matsubara frequencies can be easily eval- 
uated and yields 



AM = 



H B m 4 U 2 2 

ATT 3 k 2 F 
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coth( 
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where 
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m 4 U 2 



T 
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We see from Eq. {7J that AM increases in a field by two 
reasons. First, the field leads to a finite magnetization 
at T = 0, and second, a finite field populates the system 
with spin waves precessing at the energy fi B H. Differen- 
tiating JJJ again over H, we obtain the spin susceptibility 
in the form 



A X (T, H) = X (T, H) - x(0, 0) = ^ 2 B ATf x {x), 
where 



sinh (x) 



[sinh(2x) 



For vanishing H, i.e., at x — > 0, f x (0) = 1, and 
A X (T, H) = X (T, 0) - X (0, 0) = n%AT. 



(9) 



(10) 



(11) 



This coincides with the earlier result 
limit of large x, f x (x) ~ 2a; , and 



In the opposite 



A X (T,H) = 2fi 2 B AT\x\ = 2^ B A\H\ 



(12) 



We see that at high fields, the spin susceptibility scales 
as \H\, i.e., is non-analytic in H. 

In Fig. 2 we plot the susceptibility as a function of 
temperature at a given H , and as a function of the mag- 
netic field at a given T. Note that at a finite H, the 



FIG. 3: A7 as a function of temperature at fixed magnetic 
field (a) and as a function of magnetic field at fixed temper- 
ature (b). In the left panel, T and A X / ' HbA are in units of 
\ibH. In the right panel /ibH and Ax/fJ-%A are in units of 
T. In these units, Aj/A = f~,(fj, B H). 



Bose term in Eq. (JSJ gives rise to a negative derivative of 
OAx/dT. This in turn gives rise to a shallow minimum 
in the temperature dependence of Ax(T, H). 

The specific heat A 7 (T, H) = 7 (T, H) - 7 (0, 0) is ob- 
tained by differentiating Eq. GJJ twice over T. At H = 0, 
Aj{T,H) = -6ATC(3) EEElIll- At a finite H, 



A 7 (T, H) = -6ATC(3) + 2AT 
{x cosh x — sinh a;) = —ATf 7 (x) 

where 

Mx) = 3 (Li 3 (e- 21 
+6C(3) - 2a; 3 + 4a; 3 cothx 



^ H ' T dxx 3 



sinh 3 x 



(13) 



2xLi 2 (e- 2x ) 
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sinh 2 x 



2x 2 log(l 
(sinh 2x - 



-2x 



)) 



and Li are polylogarithmic functions. At a; <C 1, fy(%) ~ 
6C(3) - ^ and 

A7(^«1)--At(6 C (3)4(^) 4 ). 

In the opposite limit of a; 3> 1, / 7 (x) = 3£(3) + 4x 4 e~ 2x , 
and 



A 7 (^ » 1) 



-AT 3C(3) + 4 



-2h b H/T 



We see that in both limits the temperature dependence 
of the specific heat is linear in T, but the prefactor 
changes by a factor of 2 between small and high fields. 
This result could be anticipated as a high magnetic field 
eliminates the non-analyticity in the polarization bubble 
IV*(q, fl n , ), such that only the second term in Eq. (J2J 
contributes to the T term in A 7 (T, H). 

The extension of the above results to an arbitrary U(q) 
is straightforward. For the susceptibility, the prefactor in 
Eq. ©, contains U{2k F ) instead of [/jHHGil- For the 
specific heat coefficient, we have, instead of <|13|) 



A 7 (T,i/) = -^Krx 



27T 3 fc! 



U(0) - -U(2k F ) 



U 2 {2k F ) 



, / 7 (x)-3C(3) \ 
■ 3C(3) J 



(14) 
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The combinations U(0) - l/2U(2k F ) and -l/2U(2k F ) 
are charge and spin components of the scattering ampli- 
tude A(tt), respectively. At large x, f-y(x) ss 3£(3), and 
the last term in the r.h.s. of l|14l) vanishes. This ob- 
viously implies that at a large field, only the charge the 
longitudinal spin components of the scattering amplitude 
contribute to A7(T). 

We next consider the 3D case. The polarization oper- 
ators at small q are given by 



n n {q,n n ,o) = - 

U^(q,n n ,T) = 



2^ 



arctan 



(Q„ — i S/i) 



v F q 



2vr 2 



arctan 
v F q v F q 



^- + ... (15) 



As before, dots stand for analytic terms, expandable in 
powers of f2 2 or q 2 , and 5fi = /ij — /ij = 2fi B H. Dif- 
ferentiating the thermodynamic potential, Eq. |0J, with 
respect to H, we obtain 



AM 
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dH 



4ir 5 v F 
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Differentiating further with respect to H, we obtain 



A Xs (T,H) = 



2n 5 v F 



M 



A I 



, (16) 



where M ~ Ep /T is the upper cutoff in the summation 
over frequency. The first term in the r.h.s of Eq. (|l(j|) is 
the susceptibility at zero field. By power counting, one 
might expect the spin susceptibility Xs (T) in 3D to scale 
as T 2 logT. However, the Matsubara sum ry^ n=1 f2„ 
only contains a T-independent term, of order E F , and a 
term — (l/6)irT 2 . This last term is universal, but it is an- 
alytic in T. As a result, Ax s (T, 0) cx T 2 is analytic and 
essentially irrelevant as the analytic in T contributions 
to Xs(T) are already present in the Lindhard function 
for free fermions. The absence of the non-analytic tem- 
perature correction to the spin susceptibility in 31? was 
first noticed in Ref. 0, (see also H). The second term in 
the r.h.s. of (|16|) is the extra contribution in a finite field. 
Evaluating the Matsubara sum we find that this contri- 
bution scales as H 2 log{max(T, /zs-ff)}. We see therefore 
that in a finite magnetic field, Xs(T) does indeed become 
non-analytic. Casting Ax s (T,H) into the scaling form, 
we obtain 



A Xs (T,H) = xo 



mUkp \ ' 
~2^ 2 ~) 



T 
Ef~ 



.(eg). (") 



where xo = ^b^f I '(^vr 5 ''E lp) is Pauli susceptibility, and 
to a logarithmic accuracy, 



g(x) = x 2 log 



E f 



T * max{x, 1} 



(18) 



The H 2 \ogH dependence of x{H) was earlier reported 
by Misawa (l6j . However, his prefactor is different from 
the one we obtained. 

Differentiating the thermodynamic potential twice 
over T, we also obtain field dependence of the specific 
heat coefficient. The field dependence in 3D parallels the 
one for 2D systems. Namely, at zero field, 



A 7 (T,0) 



3 [mkpY 

20 7T 2 



(*7(0) - -U{2k F )) 2 + -U 2 {2k F ) 



T 
Ep~ 



In 



(19) 

Ep_ 
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In a finite field, the charge part is not affected, while in 
the spin part, the logarithmic factor 3 log ^f- is replaced 



bylog^ + 21og- 



T m ax{T,„ B H\ - As a reSult > at VE>H > T, 

Aj(T, H) still behaves as T 2 logT, but the prefactor gets 
smaller. 

To summarize, in this paper we analyzed non-analytic 
terms in the magnetization, the spin susceptibility and 
the specific heat of 2D and 3D Fermi liquids, placed into 
an external magnetic field /is-ff -C Ep. We obtained the 
non-analytic terms in the forms of scaling functions of 
HbH/T. We found that at hbH 3> T, the spin suscep- 
tibility scales as \H\ in 2D and as H 2 log \H\ in 3D. The 
specific heat in a field preserves the same temperature 
dependence as in the absence of a field, but the prefactor 
changes between small and large fisH/T. 
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